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Abstract 

The spin-glass g-state Potts model on d-dimensional diamond hierarchical lattices is investigated by an exact real space 
renormalization group scheme. Above a critical dimension di(q) for q > 2, the coupling constants probability distribution 
flows to a low temperature strange attractor or to the high temperature paramagnetic fixed point, according to the 
temperature is below or above the critical temperature T c (q,d). The strange attractor was investigated considering 
four initial different distributions for q = 3 and d = 5 presenting strong robustness in shape and temperature interval 
suggesting a condensed phase with algebraic decay. 
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1. Introduction 

The g-states Potts model, proposed a longtime ago by 
Domb as the subject of Potts doctoral thesis [l[ , has found 
a wide range of applicability in many fields of both basic 
and material sciences. The Potts model was conceived as 
a generalization of the Ising model |2J, when q = 2, and 
the Askin- Teller model {q = 4) [3J. It also mimics the 
problem of percolation (q = 1) [4j, [5j and even the prob- 
lem of the linear resistor networks (q = 0) |6j . All of the 
above mentioned problems were also encompassed by the 
random-cluster model introduced by Fortuin and Kaste- 
leyn [Tj. Furthermore, the degeneracy of the ground state 
of the antiferromagnetic Potts models was shown to be re- 
lated with the q- coloring problem [7|. It is important to 
emphasize that the most important feature of the math- 
ematical structure of the Potts model is the equivalence 
between its partition function and Tutte polynomial [8|. 
Concerning applications, the Potts model has been ap- 
plied in many fields, such as biology |9j , sociology [10J and 



material science [llj. In the latter, for instance, the tech- 



nique of Monte Carlo simulations on the Potts model has 
been applied to a wide variety of phenomena, such as diffu- 
sion in polycrystalline microstructures 12| and the study 



of viscous instabilities in foam-flow behavior [131 ] . 

In this letter, the properties of the g-state Potts model 
with random competing interactions are investigated. This 
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model is called Potts glass in allusion to the particular case 
when q = 2, widely known in the literature as the Ising 
spin-glass model. The absence of spin-inversion symmetry 
and a different nature of the frustration effects distinguish 
it from its Ising counterpart, exhibiting rather a richer crit- 

M,E^- The pure 



ical behavior in mean-field theory [14j, 
and disordered Potts model (without frustration) and re- 
lated models have been widely investigated in the past. 
For review papers see, for instance, [6j, [17J. The present 
work focuses on the study of the g-state Potts model with 
random frustrated exchange interactions on a family of 
diamond-type hierarchical lattices [18( with scale factor 
6 = 2. When symmetrical zero-centered random exchange 
couplings distributions are considered the system under- 
goes a phase transition from a paramagnetic high temper- 
ature phase to a low-temperature condensed phase above 
some dimension, di(q). The phase diagram of the Potts 
glass model had been investigated before, by the Migdal- 
Kadanoff (M-K) renormalization group (RG) scheme, in- 
dicating the presence of a condensed phase at finite tem- 
peratures when q > 2 and d = 4 19]. In reference |19|. 
however, it was assumed (as working hypothesis) that the 
initial symmetrical Gaussian probability distribution of 
coupling constants with variance a is transformed under 
renormalization into another symmetrical Gaussian distri- 
bution but with variance \/b d ~ 1 tj. 

The M-K real-space RG method for Bravais lattices 
is known to be equivalent to exactly solving the model 
on diamond-like hierarchical lattices 2JJ, [21 1 . In a recent 
paper [22j, however, the authors studied the ocurrence 
of phase transitions of the Potts glass model using this 
exact approach by numerically following the flow of the 
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renormalized probability distribution in an appropriated 
parameter space. Such space was previously considered in 
reference [23J, [2J] to study the Ising spin-glass model. The 
q-state Potts glass model was considered on lattices with 
several fractal dimension, determining the critical temper- 
ature and the upper and lower bounds for the associated 
lower critical dimension di(q). For instance, for q = 3 
the lower (upper) bound was found to be 4.46 (4.58), in 
contrast with the result obtained in reference[19], which 
founds the transition occurring for d < 4. Here we fur- 
ther explore the flow of the renormalized probability dis- 
tribution in the whole parameter space and investigate the 
nature of the low temperature stable fixed point, which 
surprisingly appeared like a strange attractor. 

2. Renormalization Procedure in Disordered Sys- 
tems 

For pure systems, the renormalization procedure con- 
sists in finding the equivalent exchange interaction for a 
pair of spins after eliminating several spins in the lat- 
tice. For a disordered system, however, the renormaliza- 
tion procedure will affect the whole distribution of cou- 
pling constants, the renormalized distribution, P'(J), is 
related with the previous (non-renormalized) distribution, 
P(J), by, 



P'(J) 



/■■■/ II P(Jio)dJiAK-K'(K)), (1) 



where K = j3J is a reduced coupling constant, K'{K) is 
the renormalization equation, and the product runs over 
all the pairs of spins < ij > . 

Equation (JTJ) should be iterated until the renormalized 
distribution reaches a fixed point distribution, character- 
istic of the thermodynamic phase. A zero-centered Dirac- 
dclta distribution, for instance, indicates a paramagnetic 
phase. The procedure adopted in this work is to produce a 
sample of random coupling constants from an initial prob- 
ability density function, feed the renormalization equation 
to find a sample of the same size of renormalized coupling 
constants, estimating the physical quantities numerically 
from the samples. The process is repeated until the fixed 
point distribution is reached. 

3. The Potts Hamiltonian and the Renormaliza- 
tion Equation 



The Potts Hamiltonian is written as, 



u 



X qJi ~ 



(2) 



<»j> 



where the sum is taken over all bonds in the lattice, the 
Kronecker 5 symbol takes the values 1 if <7, = ctj or 
otherwise, and <7j = 1, 2, . . . , q are the g-states Potts spins 
variables, located at each site of a diamond-type hierarchi- 
cal lattice with p branches and scale factor b. The basic 
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Figure 1: Renormalization group scheme on the dj dimension 



diamond-type hierarchical lattice with p 
tor 6 = 2. 



/ 
branches and scaling fac- 



unit of such lattice is illustrated in figure [TJ where \i and \j! 
are called external sites and the set {<Ji} represents the in- 
ternal sites 18j. Ki and Li are reduced coupling constants, 



Ki = PJ^i and L, = fij^i^. The lattice generations or 
hierarchies are successively built by replacing each connec- 
tion of the basic unit by the basic unit itself, yielding to a 
graph with fractal dimension, 



d f = l 



lap 
ln~2' 



The exact renormalization process on a n-generation 
lattice consists in partially tracing the partition function 
along all the internal sites introduced in the n th generation 
leading to a (n-1)— generation lattice with a set of effective 
reduced coupling constants {K^} given by 



1 p 



(q- l)+cxp(qKi + q Li) 



(q — 2) + exp (q Ki) + exp (q Li) 



(3) 



Equation ([3]) is the local renormalization equation. The 
right-hand side of equation ([3]) receives values for the re- 
duced coupling constants calculated from the n-generation 
of the coupling constant distribution, P(J), resulting in 
one of the possible values of the renormalized reduced 
coupling constant of the (n — regeneration. The renor- 
malization procedure starts from the thermodynamic limit 
(generation n — > oo) where the coupling constants are as- 
sumed to have a well known distribution, actually, Gaus- 
sian, Delta-Bimodal, Uniform, or Exponential, and well 
defined temperature T. 

3.1. Probability distribution renormalization flow 

For each n-generation lattice we can define a set of 
thermal transmissivities variables {Uj}, each one associ- 
ated with the respective bond, ie., 



U 



1 - cxp(-qPJij) 
l + (q- l)cxp(-qf3Jij)' 



(1) 



Thermal transmissivity tij represents the pair correlation 
function I\j between sites ij |17| . 

A system with a probability distribution of coupling 
constants P(J) = P ({Jij}) yields a thermal transmissivity 



variance, 

A 2 = [fa - [U 3 ]jf]j 



(5) 



where [• • •] j means the average over the probability distri- 
bution, P(J). 

Every fixed probability distribution has a signature in 
the diagram A 2 x T when the temperature is varied. In 
the case of the initial probability distributions considered 
in this work such representation can be seen in figure [5] 
for the case q = 3. Notice that the T — ► limit is the 
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Figure 2: A 2 behavior for distinct probabilities distributions against 
temperature for the q = 3 case. 



same for all cases indicating an universal zero tempera- 
ture behavior. Such universal point (A* 2 ,T = 0) can be 
analytically obtained in the case of the Bimodal distribu- 
tion, P(Jy) = [S(J l3 - 1) + 5(Jij + I)]/2, ie. 



A 



*2 



2(9-1) 



(0) 



Note that for q = 2, the Ising limit case, A* — 5- 1 is recov- 
ered as seen in [23j. 

For each application of the renormalization equation, 
an effective temperature proportional to the inverse of the 
square root of the reduced coupling variance may be de- 
fined, ie. 



T r ex 



1 



y/KKij - [K^J^J 



(7) 



In this work, the flow of the renormalized probabil- 
ity distributions is followed numerically in the parameter 
space A 2 x T r , each point representing a stage of the renor- 
malization process. 



3.2. Strange attractor 

Four distinct initial symmetrical probability distribu- 
tion were considered, namely the Gaussian, the delta-bimodal, 
the uniform and the exponential ones. For lattices with 
dimension greater than di {q) and independently of the na- 
ture of the initial distribution, the renormalization flow 
displays the following same features: the paramagnetic 
phase is characterized by a fixed point at infinity renormal- 
ized temperature and zero transmissivity variance, while 
the low-temperature phase is characterized by a strange 
attractor in the renormalization flux for q > 3 - or a fixed- 
point A = 1, T r = for q = 2 as previously reported [23j 
for the Ising SG. The observation of such strange attractor, 
which is located in a region of low but finite temperatures, 
is the main result to be reported in the present work. 

A first signature of the existence of this strange attrac- 
tor was reported by Banavar and Bray [25] long time ago. 
These authors observed an unusual behavior in the renor- 
malization of [Kij] j and its standard deviation in the case 
d = 5 and scale factor 6 = 2 (present case). In their words 
these quantities "wanders chaotically in a smallish region 
around [Kij]j ~ 14 ... (persisting) for 1000's of iterations 
and happens for a wide range of starting temperatures" . 
Furthermore, they do not give a clear explanation of this 
behavior but speculate that the delicate balance between 
energy and entropy causes that wandering phase in d = 5. 

This study uses the same methodology used in refer- 
ences [23j,|25(, ie. the successive construction of large pools 
of renormalized constant couplings, calculated using equa- 
tion ([3]) starting from an initial probability distribution 
and fixed temperature. Figure [3] illustrates the renormal- 
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Figure 3: Renormalization flow in the A 2 X T parameter space for 
the q = 3-Potts model in d = 5 dimension for an initial symmetric 
Gaussian distribution. Dotted curve indicates the low-temperature 
flow while continuous one illustrates the high-tempcrature flow, ac- 
cording to the arrows. Inset magnifies the flow bifurcation region. 

ization flow diagram in the A 2 x T space for the q = 3 
case in d — 5 dimensions, starting from a Gaussian distri- 
bution of the exchange couplings in two temperatures. For 
initial temperatures below ~ 3.0 the renormalized distri- 
bution of interactions (solid triangles) flows to the strange 



attractor located in a finite region of the parameter space, 
for high initial temperatures (T ~ 3.5), the renormalized 
distribution of interactions (solid circles) flows to high tem- 
peratures and zero thermal transmissivity variance, char- 
acterizing the paramagnetic phase. Figure S] displays the 
strange attractor constructed with 10 5 steps of renormal- 
ization, lying in the interval (0.15, 0.35) starting from two 
distinct initial distributions, (a) the Delta-Bimodal and 
(b) the Gaussian ones. In reference [22|. some of the present 
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Figure 4: Attractor for two distinct initial distribution, (a) Bimodal 
and (b) Gaussian 



authors, using the same methodology, found the critical 
temperatures of the model for several lattices, indicating 
the ones where the transition may occurs for appropriated 
values of d and q. 

The fractal dimension of the attractor was estimated 
by box counting, yielding to, 

D= 1.66(1), 

for every initial distribution and for the case q = 3 and 
d = 5. A proper analysis of the attractor is taking place 
at the moment, by estimating its Lyapunov exponents by 
the Wolf method [26[ with various values of q and d. 

The strange attractor describing the condensed phase 
occurs at a non-zero temperature interval probably due 
to the existence of residual entropy, such as the unusual 



algebraic phase previously studied [27|, [28J, l29| for the anti- 
ferromagnetic Potts model. Detailed studies of the correla- 
tion functions and the local magnetization are under devel- 
opment in order to determine the nature of such condensed 
phase. Preliminary studies of the present model on lattices 
with scale factor 6 = 3, which corresponds to bipartite 
lattices, indicates however the absence of the strange at- 
tractor in any dimension. Moreover, the condensed phase 
which occurs above a certain lower critical dimension di (q) 
does not presents characteristics of spin-glass phase |30jj - 

4. Conclusions 

The q-state spin-glass Potts model undergoes a phase 
transition whenever d > di(q), from a high temperature 
paramagnetic phase to a low temperature condensed phase, 
characterized by strange attractors. These strange attrac- 
tors are superimposed in the same region of parameter 
space independently of the initial distribution of coupling 
constants and initial temperatures below T c , even if this 
initial temperature is chosen close to zero, below the at- 
tractor region. Such strange attractors, however, change 
in shape and position for distinct {q, d) models in lattices 
with scale factor b = 2. The occurrence of the strange at- 
tractor in a finite non-zero temperature interval suggests 
a condensed phase with algebraic decay of the correlation 
function. 
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